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Atom interferometer as a selective sensor of rotation or gravity 
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In the presence ol Earth gravity and gravity-gradient forces, centrifugal and Coriolis forces caused 
by the Earth rotation, the phase of the time-domain atom interferometers is calculated with accu- 
racy up to the terms proportional to the fourth degree of the time separation between pulses. We 
considered double-loop atom interferometers and found appropriate condition to eliminate their sen- 
sitivity to acceleration to get atomic gyroscope, or to eliminate the sensitivity to rotation to increase 
accuracy of the atomic gravimeter. Consequent use of these interferometers allows one to measure 
oXj, all components of the acceleration and rotation frequency projection on the plane perpendicular to 

3 ' gravity acceleration. Atom interference on the Raman transition driving by noncounterpropagating 

■^L I optical fields is proposed to exclude stimulated echo processes which can affect the accuracy of the 

atomic gyroscopes. Using noncounterpropagating optical fields allows one to get a new type of 
the Ramsey fringes arising in the unidirectional Raman pulses and therefore centered at the two- 
quantum line center. Density matrix in the Wigner representation is used to perform calculations. 
It is shown that in the time between pulses, in the noninertial frame, for atoms with fully quantized 
r^ ' spatial degrees of freedom, this density matrix obeys classical Liouville equations. 

I. INTRODUCTION 

d Since atom interference [ij has been proposed as a sensor of inertial effects ;2| and the use of Raman transition 

c/3 , between atomic hyperfine sublevels Q allowed tremendously increase the time separation between optical fields, un- 
O ' precedented accuracy in the measurement of the Earth rotation [Jj , gravity gradients i3 j ft/'nics iS] j ^nd accelerations 
c/2 i [3 has also been achieved. The theoretical analysis [a. l9l llOl lll| showed that the current level of interferometer phase 
^~^' measurements is sufficient to sense each source changing the atomic motion, i.e., rotation, acceleration, acceleration 
1-^ ^ gradient, and recoil effect fl^ . as well as the interplay between them, such as between rotation and acceleration [lO| . 
I '~^ . recoil effect and rotation [l3|, recoil effect and gravity gradient ]1(|| . 

At the same time for navigation and geodetic applications one needs to measure separately a gravity acceleration g 

^^ , and rotation frequency Jl, i.e., one needs an interferometer whose phase is selectively sensitive only to one component 

^ ' of these vectors, while the sensitivity to others can be either excluded with sufficiently high accuracy or precisely taken 

^^ , into account. An example here is an atomic gyro j^], where the Jl projection on the gravity acceleration g has been 

measured, because influence of gravity has been excluded using spatially separated fields propagating in horizontal 

plane and using signals from two counterpropagating atomic beams. 

f~^ . In this article we consider the theory of atom interferometers that would allow one to measure separately an atom 

\^ ' acceleration g and the rotation frequency component il which are perpendicular to g. We explore the fact that, 

^^ , when rotation, gravity gradient and recoil effect only slightly affect the atoms' trajectory, parts of the interferometer 

(/2 phase associated with g and il evolve correspondingly as T^ and T^ while other contributions are precisely known or 

O : negligibly small. Using particular kinds of the atom interferometers one can exclude one of these dependences getting 

c/j a sensor of acceleration or rotation only. 

J>~^ Usual time domain atom interferometer, consisting of three resonant pulses applied to the cold atom cloud at 

1-5 moments t = 0, T2 and T3 could not serve for these purposes because the only parameter, ratio of the time separation 

between pulses T3/T2, is already used to eliminate linear in time part of the phase. This Doppler phase vanishes at 
the echo point, Tk = 2T2. To eliminate another terms one needs at least a four- pulse interferometer. It was found 
previously y, llJ, uM ^^at double loop interferometer consisting of two 7r/2 pulses at the beginning and end and two 
TT pulses in between, i.e., the (7r/2,Ti = 0) — (tt, r2) — (7r,T3) — (7r/2,T4) interferometer, has no sensitivity to the 

. -1 o 

. . 1 homogeneous acceleration if T2 — —T4 and T3 — -T4, see fig. H^a). 

One can expect that the double- loop interferometer can serve as a sensor of rotation. It is correct only if pulses' areas 
are precisely equal to assigned values. If one could not hold areas equal to tt or tt/2, then stimulated echo processes 
also contribute [see Fig. Wi^)] to the interferometer phase and since stimulated echo is insensitive to the atomic 
motion between second and third pulses, their contribution is acceleration sensitive. Since typically contribution to 
the interferometers phase caused by the acceleration is 5-6 orders of magnitude larger than that caused by rotation 
[iJI, it is extremely important to find a technique for excluding stimulated echo. To achieve this goal one should 
violate the condition for stimulated echo. This opportunity exists if one uses Raman transitions, where the splitting 
of atomic momentum states occurs under the action of two optical fields having wave vectors qi and q2, such as 
an effective wave vector associated with Raman transition is given by k = qi — q2. In principle one can use for this 
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FIG. 1: Recoil diagrams for four-pulse atom interferometers, (a) double- loop interferometer insensitive to homogeneous ac- 
celeration; (6) stimulated echo processes, shown by black parallelograms, which could restore sensitivity to acceleration; (c) 

double- loop interferometer with nonidentical loops, T2,s = T4. Solid and dashed lines correspond to the lower and upper 



atomic states. 



urpose noncouterpropagating fields. Interferometers witli qi 7^ — q2 liave been previously created jig and studied 
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Since for Raman transition between hyperfine sublevels Qi ~ 92, |k| ~ A: sin a, where a is tlie half-angle between 
vectors qi and q2 , and k w 2qi is the maximum value of the effective wave vector corresponding to counterpropagating 
optical fields. If kj is an effective wave vector associated with Raman pulse j, than the phase matching condition for 
pulse sequences 7r/2 — it — it — 7r/2 [Fig. CJa)] and 7r/2 — 7r/2 — 7r/2 — 7r/2 [Fig. n[b)] are 



k4 - 2k3 + 2k2 - ki == 0, 

k4 - kg - k2 + ki == 0. 



(la) 
(lb) 



These conditions mean the same only for parallel pulses having effective wave vectors of the same absolute values kj . 
But even for pulses with slightly different kj one can hold only one of this condition and violate another. In this paper 
we propose and analyze an exclusion of the stimulated echo processes using for double loop atom interferometers 
Raman pulses with nonequal effective wave vectors. 

While eliminating T^ dependences one can use time domain interferometer as a gyroscope, to measure rotation 
frequency components perpendicular to acceleration, eliminating T^ dependences would guarantee that contributions 
of this order, which can be caused by J10j | gravity gradient (for atoms with nonzero launch momentum), rotation 
(when effective wave vectors are not parallel to the acceleration), recoil effect and combinations of these factors, do 
not affect the interferometer phase at all. Absence of T^ dependences should increase the accuracy of the acceleration 
measurement. We will show below that T^ dependences disappear for the double- loop interferometer with nonidentical 



loops shown in Fig. ^c). Since in this case pulse separations T2, Tj, and T4 are incommensurate, one can be sure that 
no other echolike processes can affect the interferometer phase. 

In this paper we perform calculations of interferometers' phases for three- and four-pulse cases, for arbitrarily 
directed gravity acceleration and rotation including Earth gravity gradient as a small correction. It would allow us to 
get the phases' terms evolving up to T** dependences. In our calculation we use equations for the density matrix in the 
Wigner representation, which are convenient to describe atomic clouds' evolution between pulses [13, because we show 
that in this representation density matrix evolves exactly as the classical distribution function in the rotating frame 
in the presence of the homogeneous gravity and gravity gradient terms. This technique is mathematically equivalent 
to the pure quantum consideration of the atomic wave function evolution in the space between pulses elaborated in 

0. 

Using density matrix in the Wigner representation allows us to evaluate precisely all effects related to the atomic 
spatial motion quantization without calculating path integrals 19] . 

The paper is arranged as follows. A simplified consideration of four-pulse interferometer is performed in Sec. ^] 
Evolution of the density matrix in the Wigner representation is considered in Sec. IIIII Atom trajectory in the Earth 
rotating frame, including gravity-gradient terms as perturbation, is calculated in Sec. IIVI Density matrix jumps under 
the pulse action are obtained in Sec. In Sec. I VII we calculate the phase of the 7r/2 — tt — 7r/2 pulse interferometer. 
Section IVIII is devoted to study four- pulse interferometers as atomic gyroscope or gravimeter, while in Sec. IVIIII we 
summarize our consideration and discuss further possible developments of the proposed multiple loop interferometers 
technique. 

II. FOUR-PULSE ATOM INTERFEROMETERS 

Consider an interaction of cold atoms with a sequence of four Raman pulses, 7r/2 — tt — tt — 7r/2, having the same 
effective wave vector k and applied at moments Ti = 0, r2, TJj, and T4. We will show below that for atoms with initial 
coordinate and momentum r and p the phase of interferometer associated with atomic motion is given by 



= k- 



Kk \ ^ f Kk \ ^ f hk 

R(r,p+— ,T4J -2Rfr,p+— ,T3J +2Rfr,p+— ,T2 



(2) 



where the function R (r, p, t) is the position of an atom at the moment t under the classical motion in rotating frame. 
For small rotation frequency ft {fit ^ 1) the atom trajectory is given by 

Rir,p,t) ^ r + ^t + g*- - n X (r+^t + gMt, (3) 

TO 2 \ m 6 J 

where to is the atomic mass, and for the phase |(5J| one arrives at the expression 

hk 




^ (k X g) . nrl (1 - 2il + 2tl) , (4) 

where 

h = T,/T,. (5) 

If the purpose is to measure the rotation frequency, then one may set to first and second terms in Eq. Q). System 

1 - 2^3 + 2t2 = 0, (6a) 

l-2tl + 24 = (6b) 

has a solution ^2 = 1/4,^3 = 3/4 corresponding to the double loop interferometer |2,llJ,ll3 shown in Fig. QJa). In 
this case 

0^ = ^ (k X g) . nri. (7) 



If the gravity is measured, then the phase ^ can be used to measure perpendicular to gravity components of 
the rotation frequency. For example, for Cs interferometer having the effective wavelength 426 nm and total time 
separation between pulses T4 = 0.8 s, at the latitude 41° the maximum value of the phase Q 0„ sa 250. 

If the purpose is to eliminate cubic terms in order to increase the accuracy of the acceleration measurement, one 
should choose the solution of the system 

1 - 2^3 + 2t2 ^ 0, (8a) 

l~2tl + 24 = 0, (8b) 



which is 



t..-^, (9) 



and for atoms having zero launch momentum, p ~ 0, the phase is given by 

0™ = -^^k.gr^ (10) 

For a given maximum separation between pulses, the absolute value of the phase IjlOII is approximately 4.2 times less 
than the signal in the usual three-pulse interferometer. Therefore, it is a kind of trade off, when one should decide 
whether it is more important to get a larger signal or to eliminate T^ corrections to the phase. 

III. DENSITY MATRIX EVOLUTION IN THE FREE SPACE 

Starting with the Lagrangian for an atom in the frame rotating with constant rate, 

2 

TfW Til o 

i:=-— + mv(flxr) + -(^xr)^-C/, (11) 

where U is the potential, one gets the generalized momentum and Hamiltonian, 

p = TOv + 77T, (Jl X r) ; (12a) 

p2 

H =-^— + 'p-{rxn) + U. (12b) 

2m 

In the space free of the laser field, i.e., free space, the density matrix in the Wigner representation |2(]| 

p(r,P,i) = 3 dsp(r+-s,r- -s, t j exp (-ip • s/?i) , 

where 

p(r,r') = i/'(r)®i/;t(r') (14) 

is the density matrix in the coordinate representation, evolves as 



(13) 



dtp (r, p, t) = / ds exp (-ip • s/h) \ 

iTi{2iTh) J {, 

2m 



2m 



P(ri,r2,t)j. . (15) 

ri=r+is,r2=r— is 



Using the fact that 



9ri = -^r +9s; 

dr, = ^dr - d, (16) 



and replacing dg => ip/Ti, one arrives at equations 

{9t + (- + r X n^ dr + [(p xn)- drU] dp + Q^p (r, p, t) 



Q = -{ihy 



U\v+ iiaSp j -U\v- i^a5p 



drUdr. 



(17a) 
(17b) 



When the potential U is not higher than bilinear function of coordinates, quantum term disappears, Q = 0. Since the 
Hamiltonian equation for the classical momentum and coordinate is 



in 



- r X $7; 



one concludes that, if 



-drU + (p X J7) ^ 



Q^n'{dlu)dl<^T-\ 



(18a) 
(18b) 

(19) 



where T is a typical time of the system's evolution, then with an accuracy QT the quantum density matrix in the 
Wigner representation obeys the same equation as the classical density matrix 

(at + far + pap)p(r,p,t) = o. (20) 

If the density matrix is known at some previous moment t' < t then the solution of Eq. 12U|) is given by 

p (r, p, t) = p [R (r, p, t' - t) , P (r, p, t' - t) , t'] , (21) 

where functions R (ro,po, t) , P (ro, po, i) arc solutions of Eqs. (|18|) subjected to the initial condition (r,p)j^Q = 
(ro,po) • Clearly these functions have to satisfy multiplication law 



R [R (r, p, ii) , P (r, p, ii) , ^2] = R (r, p, ti + ^2) : 
P [R (r, p, ii) , P (r, p, ii) , t2] = P (r, p, ii + ^2) • 



(22a) 
(22b) 



Solution (|21() is the consequence of the phase-space invariance in time. This is correct only for canonical variables 
(r, p) and incorrect for variables (r, v) . Solution H21|l can be proved as follows 



(9p(r,p,t) _ dp dRi{r,p,t' -t) dp dPi{r,p,t' - t) 



Since 



and 



one gets 



dt dR, 




dt dP^ 


dt 


dR, (r, p, t' - 


-t) 


dRi{r,p,t' -t) 


dH 


dt 




dt' 


dp 


dPi (r, p, t' - 


-t) 


dP,(r,p,t'-f) 


dH 


dt 




dt' 


dR^' 



dp (r, p, t) 
dt 



{H.p}r,p, 



(23) 



(24) 



(25) 



(26) 



where {H,p} is a Poisson bracket, and using invariance of the Poisson brackets {H,p}^p = {H,p} one arrives at 
Eq. (HOI). 

IV. ATOM CLASSICAL MOTION IN THE ROTATING FRAME 



In this section we calculate the atomic classical trajectory in the potential consisting of the linear part and gravity 
gradient part 



U (r) = —ma. ■ r + SU (r) , 



SU{r) = ^T,, 



3{n-ry 



(27a) 
(27b) 



where a = g — J7^i?^'X {u x n) is a sum of the gravity and centrifugal accelerations, v — $7/17 is a unit vector along the 
rotation frequency, n = R/i? is a unit vector along the system's displacement R from the Earth center, T^^ = ~'^9z/ R, 
Qz is the vertical component of the gravity acceleration. 

We will treat the gradient term Ij27bp as a perturbation. In the zero order the atom evolves as 



f = r X $7 H , 

ra 



p= (p X i7) + ma, 

subject to initial conditions r (0) = ro, p (0) — po, while first-order corrections are given by 

m 
5p = -^TzzV^, 



Sr = --TzzY^, 



where r™ and p^ evolve as 



r-y = r-y X J7 + p^ , 

P7 = P7 X ^ + r — 3n (n • r) , 

subject to initial conditions r^ (0) = 0, p^ (0) — 0. 

One sees that to get the atom trajectory it is necessary to solve equations of the type 

b = f^b X 1/ + f (t) . 

Directing temporarily z axis along v, one finds that variable S, = b^ — iby is given by 



^ = exp {iQt) £,Q + dr exp [iQ 
Jo 



it~T)][fAr)-^fyiT)], 



so that the solution of Eq. (|31|l can be found as 

b = cosnt{bQx,boy,0) +smnt{bQy,~bo^,0) 



+ / dT {cos n{t~T) [U (r) , fy (r) , 0] + sin f! (f - r) [fy (r) , -/, (r) , 0]} + 



(28a) 
(28b) 

(29a) 
(29b) 



(30a) 
(30b) 



(31) 



(32) 



0,0,60.+ / drfzir) 



(33) 



Since for the chosen coordinate system {box,boy,0) = —ux {v x b) and (6oy, — feoa;, 0) = — z/xb, then in the vector 
representation the solution of Eq. H31|) is given by 



b = b/i + bp, 

b/i = f X (i/ X bo cos Q,t + bo sin Vtt) + v {v ■ bo) 



(34a) 
(34b) 



(34c) 



(35) 



hp= drj-i/ X [1/ X f (t) cosr2(i- t) + f (t) sin f^(t-r)] +iy[i/-f(T)]}. 
Jo 

If in turn the driving term in Eq. 131() can be represented as 

f{T) = vx[vx sqi (r) + sq2 (r)] + u {u ■ s) qo (r) , 
then 

hp = V X <u X / drs [qi (r) cos J7 (f — r) — 52 (t) sin i7 (< — t)] 

+ / drs [qi (r) sin fi (t — t) + 52 (t) cos fl(t — t)]> +u dr {v> ■ s) qo (r) . 

Consecutive application of Eqs. H34|l - H36() brought us to the following expression for the function R(r,p,t) 

R (r, p, t) = Ro (r, p, t) + 611 (r, p, t) ; (37a 



(36) 



Ro (r, p, t) = cos X (r ^ 1] — sin xu x (r -\ 1] +(1 — cos x) v 

\ ra J \ m / 

/o,i (x) = x^ (cosx + xsinx — 1) , 

/o,2 (x) — x^ (sin a; — x cos x) ; 



r+ — t 

TO 



(37b) 



(5R(r,p,t) = --T^^r^; r^ == ri + . . . rs; 



(37c) 



Y-i = t^ \u u ■ { — I 1] (1 — COS x) + { — \ 1 1 cos x — V X { — I 1] sin x \ 



t" ^ a/i,i (x) + i^{a-u) 



24 



- /i,i (2^) 



- ly X a/i,2 (x) > , 



/i,i (x) = x"'^ [l - cosx (1 - 2-^x'^) - xsinx (l - G^^x^)] 
/i,2 (x) = x^ [xcosx (1 — 6^ X j — sinx (1 — 2 



-^x^l 



r2 ~ 3t (n- v>) < V (n- u) 



('^ • r) hi (x) 



1 






n 



[/22(x)-6-i]+<2i..a(^/23(x)-l 



/21 (x 

/22 (x 

/23(a; 

f2b{x 
f2&{x 



X '[3(1 — COS x) — X sin x — 2 
= x^ [2 (1 — cos x) — X sin x] , 
= x^'' (x (2 + cos x) — 3 sin x) ; 



-1x2] 



(37d) 



i^ ■ (r/21 {x) + ^t/22 (x) + atV23 {x))] +uxn\u- (r/02 (x) + ^i/25 (x) + at2/26 {x))] ] , (37e) 

\ TO /JL\ TO /J J 

x^ (cosx + xsinx — 1) , 

x^ [2sinx — X (1 + cosx)] , 
4 



rg = -4-13^2 [(ly X n) • (iv X r)] [/30 (x) u {n ■ u) + /31 (x) n + /32 (x) i. x n] , 
/30 (x) = x^2 [4 — (4 + x^) cosx — xsinx] , 
/31 (x) = x^ (sin X + X cos x) , 
/32 (x) = x^2 [x cos X — sin x (l + x^)] ; 



(37f) 



r4 



4to, 



t^ [{u xn) ■ {v X p)] 



/40 (x) zy (n • iv) 



/41 (x) n — - sinxi/ x n 



/40 (x) = X [7 sin X — x(3 + 3 x) cos x — 4x] , 
/41 (x) = x^ [sinx — X (1 — 3^ x ) cosx] ; 



(37g) 



rg = -4 ^3^2 [(1/ X n) • r] [/50 (x) v> (n ■ v>) + /51 (x) n - /02 (x) iv x n] 
/50 (x) = x^ [x (4 — cos x) — sin x (3 + x j] , 
/51 {x) — x^2 [sinx (1 — x^) — xcosx] ; 



(37h) 



rg = -^*^ K*^ X n) • p] [/eo (x) z/ (n • iv) + 3 ^ sinxn + /62 (x) iv x n] 
/go (x) = x"^ [6 (1 - cos x) - X sin X (4 + 3^-^x2)] , 
/62 (x) = x^ [xcosx (1 + 3^ X ) — sinx] ; 



(37i) 



8 

3 

r7 = --t* [{u X n) • a] [/70 {x) z^ (n • i^) + /71 (x) n + /72 {x) u x n] (37j) 



/70 (x) = a;^^ {2x [4 + (1 + Q-^x^) cos x] - sin a; (10 + x^) } 
/71 (x) = a;^^ [2a; cos X (l — 6~^x^) — sin a; (2 — a;^)] , 
/72 (x) = x^ [a; cos a; — sin a; (l — 3^ x )] ; 



rg = -^t" [(«y X n) • (iv X a)] [fs,o (x) iv (n • «^) + /g,! (x) n + /g.a (x) 1^ x n] (37k) 



where x = fit. 



/so (a;) ^x^-^jS [1- (1 + 8^^x2) cos x] - xsinx (l + 3"^x2) - 2x^} 
/s.i (x) = x^ [4 — cosx (4 — x^) — xsinx (3 — 3^"'^x^)] , 
fs.2 (x) = x^* [sinx (2 — x^) — xcosx (2 — 3~"^x^)] ; 



DENSITY MATRIX EVOLUTION INSIDE RAMAN PULSE 



Consider atoms interacting with a field of two traveling waves, 

E ^- Re |Eie*(ii-""""i*+'''i) + E,2e'^'i^-''-'^^*+'^^'^\ , (38) 

where Ej, Wj, q^, and ■0. are amplitudes, frequencies,, wave vectors, and phases of waves, respectively. We assume 
that the pulse of the field (|38() start and end times are ig and to + r, where r is the pulse duration, fields Ei and 
E2 are resonant to adjacent transitions \g) — > |z) and |e) — *■ \i) correspondingly, and g and e denote lower and upper 
hyperfine sublevels of the atomic ground state manyfold, while i corresponds to the intermediate level in the Raman 
two-quantum transition. Regarding the pulse duration we assume that being larger than inverse fields' detunings 

T»|Ar\ (39) 

where 

A ^ LOi - UJig K UJ2 - t^ie, (40) 

LUig and LUie are transition frequencies, it is also sufficiently small in respect to all other relevant time intervals 

T « min {t„ \S12r' , {kAp/m)'^ , (^feVm)"' , |k • gf^/' , {nf^} , (41) 

where J 7^ 1, 

k = qi - q2 (42) 

is the effective wave vector, 

S12 = Wi - W2 - l^eg (43) 

is the Raman detuning, Ap is the width of the atoms' distribution over momenta. 

Evolution of atomic levels amplitudes under assumption (I39|) and in the absence of acceleration and rotation has 
been considered in review |21|. From Eqs. (51) and (52) in 21], using conditions (|41|l one finds that in the Shrodinger 
representation amplitudes of levels after the pulse action are given by 



c (e, p + hk, to + r) = exp {-ir {Q^'^ + flf^) /2} i cos f - j c (e, p + ^k, to) 



-i exp {-i [</« + Si2to]} sin ( - ) c (5, p, to) } , (44a) 

c (<7, p, to + r) = exp {-zr {ilf^ + nf^) /2} \ -i exp {i [</> + (5i2to]} sin ( ^ ) c (e, p + ^k, to) 



'Cos( -)c(g,p,to)J>, (44b) 



where Q^^ = ^^1^ /4A and Q^^ = |f7g|^ /4A are ac-Stark shifts of levels |e) and \g) , f7e = - (i |d • E2I e) /h, ^g 



{i\d--E,\g)/h, 



(45) 
(46) 



= V2- "01, 

is the pulse area, ^eff — |f^e^g/2A| is the effective Rabi frequency associated with the Raman transition |e) — > \g) . 
In the derivation of Eqs. 144() we assumed that Raman detuning is much smaller than the effective Rabi frequency. 
1^12 1 <C ^eff- For the optimal pulse areas 6* ~ 1, this condition follows from the inequality 14111 . 
For the density matrix in the Wigner representation H13|) . which can be expressed through levels' amplitudes as 



Pa/3(r,P,t) = / 



— ^c (a, p + -hq, tj c* (p, p - -hq, t ] e'^J ^ 



one finds 



Pee (r, P, ta + r) = cos^ (-] p^^ (r, p, to) + sin^ (-] p„„ (r, p - ^ik, to) 



hk 



Re< isin6'exp[-i(k-r-(5i2io -</>)] Pen r,p- -z-ih 



Pgg {r,p,ta + t) = sin ( - 1 p^Jr, p + ^k, to) + cos ( - 1 p^^ (r,p, to) 

-Re<^ isin6'exp[-i(k-r-(5i2to - (l>)] Pe„ (r,P+ ^^'^o 



p (r, p, to + t) = -i sin 9 exp [i (k • r - ^i2to - 0)] 



hk 

Pee 1 r,p + — ,to 



P99 ir,p-— ,to 



+ cos^ [2) ^"3 '-'"' P' ^°^ ~^ ^^"^^ ( 2 ) '^'^^ f^* (k • r - (5i2to - </>)] p^^ (r, P, ^o) 
In the specific cases of tt and tt/2 pulses Eqs. H48|l reduce to 

Pee (i": P. ^0 + t) = Pgg (r, p - tik, to) ; 

Pgg (i": P. h + t) = p^^ (r, p + ^k, to) ; 

Peg (i", P, to + t) = exp [2i (k • r - Si2to - </))] Pge (r, p, to) , 
for 9 = TT, and 

Pee (!■> P, *0 + t) = - [Pee (l", P, ^o) + Pgg (r, P - ?ik, to)] 



Pgg (r, P, to + ^) = 2 t^ee *^'"' P + ^^' ^°) + ^59 '^''' P' ^°)] 



Tik 
Re<j iexp[-i(k-r-(5i2to- 0)] Pe„ ( r, p - — -,to 



hk 
Re<{ zexp[-i(k-r-(5i2to- </»)] p^g ( r,P + ^T'^o 



;ik 

Pee 1 I",P+— ,to 



Pen (r, P, to + t) = - exp [i (k • r - S^to - 0)] 



1 



?ik 
^99 \^^P~ -^^to 



for 9^7r/2. 



+ o {Peg (r, P, to) + exp [2i (k • r - ^lato - (/))] Pg^ (r, p, to)} 



VI. THREE PULSE ATOM INTERFEROMETER 



(47) 



(48a) 



(48b) 



(48c) 

(49a) 
(49b) 
(49c) 



(50a) 



(50b) 



(50c) 



Consider first a tt/2 — tt — tt/2 sequence of pulses applied at moments Ti =0, T2 — T3/2, T3. If initially all atoms 
are in the lower state, where their distribution in the Wigner representation is / (r, p) , 



p(r,P,0) 





/(r,p) /' 



(51) 
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then applying consequently Eqs. (|21() . H49() . H50|l one finds that after the third pulse action, the upper state distribution 
is given by 

p,, (r, p, Tg + r) = 1 {/ [R (R (r, p, T2 - Tg) , P (r, p, T^ - T3) - ^k, -T^) , 

P (R (r, p, T2 - T3) , P (r, p, T2 ^ Tg) - ?ik, -T^)] 

+ / [R (R (r, p - nk, T2 - Ta) , P (r, p - Kk, T2 - T3) + hk, -T2) , 

P (R (r, p - hk, T2 - Tg) , P (r, p - Tik, T^ - T3) + hk, -T^) - hk]} 



kr-2kR r,p 



hk 



..T.-T. 



k-R r,p 



hk 



-T, 



-<5i2 (Tg - 2T2) - (/.3 + 2(/.2 - (/.J / 



ak \ „ / hk \ hk 

R(r,p--,-T3 ,P r,p--,-T3 -- 



(52) 



This expression can be used to calculate any response associated with atoms on the upper level. It can be used to 
study the spatial shift and deformation of atomic clouds. In this paper we use it to get the total probability of atoms 
excitation 



w 



/ drdvPee (i", P, 73 + r) . 



(53) 



Choosing for each term in Eq. (|52|l arguments of the initial distribution / as integration variables and taking into 
account that phase space stays invariant under both the free atom motion and recoil kicks of the momentum, one 
arrives at the following probability of excitation: 



w 



1 

2 


1 - / drdp cos (0„ 


- <t>R- <t>Z- <t>l+'^4>2) f i^^V) 


k- 


r^ / hk \ 
R r,p + — ,r3 

. V 27 


^ f Kk \ ] 
-2R r,p + — ,r2 +r 

V 2 y _ 


■> 



<^12 (J3 - 2T2 



(54a) 

(54b) 
(54c) 



where 4>j is a phase H45(l associated with the pulse j, and (f)^ is the phase caused by atomic motion and recoil effect, 
which reduces to the Doppler and recoil phase shifts in the absence of rotation and acceleration, and (/)^ is Ramsey 
phase [2^ . 

Time independent and linear in Tj parts of the 0„ cancel owing to the coincidence of wave vectors of pulses (phase- 
matching conditions) and proper relations between pulses' separations. Other terms can be mapped out as those 
depending on r and r independent. We assume that the atomic cloud size is sufficiently small to neglect r-dependent 
terms. We accept the same assumption for nonlinear in Tj parts of the 0„ regarding the momentum distribution, 
except that we allow atoms to be launched into some initial momentum p. Thereupon one can drop integration in 
Eq. H54a|l over phase space, and (/>„ becomes the phase of the interferometer. 

All phase terms up to the order T3 are presented in TableQ] They can be obtained by multiplying by k corresponding 
vectors in Eqs. (|37|l . expanding functions fij {x) and replacing factor i" {n > 1) as 



f 



r«(i 



') 



(55) 



We would like to mention that the fourth term in this table is 1.5 times smaller than the corresponding part of the 
phase calculated in 9], which has to be corrected. 

As an example, and to verify our calculations we present in columns 2 and 3 phase terms for the case considered 
in article [lO|, when atoms are launched from the coordinate system origin and gravity acceleration, launch velocity, 
and wave vectors are all vertical. 



k=A;n, 

p = rnVzH. 



(56) 



The choice of the coordinate axes are also the same, z, vertical; y, south-north; x, west-east directions. 

If Raman pulses are comprised from noncounterpropagating fields then effective wave vectors still have to be parallel 
to one another, and for any ratio T2/T3 one can find appropriate magnitudes of absolute values of effective wave vectors 
and for T2 7^ 23/2 one gets nonzero Ramsey phase (|54cf) . 
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TABLE I: Three-pulse interferometer phase shifts. Column 1, for an arbitrary acceleration a, rotation frequency ft, and wave 
vector k; column 2, for the vertical gravity acceleration (g =gz'n) , launch velocity (p = mv^n) and wave vector (k = fcn). 
Those dependencies which appear in the different cells of column 2 are underlined and summed up in the column 3. 

jkg.Ti 

hnlRTi 



-hvzQlTi 






4 ■^ 






oo2 


^;^(«xk). 


pTi 




8m 


[{kiy) (pi/) - 


(kp)] 


Ta 




in- (kx a 


)Ti 




16m ^ 


k^] Q' 


Ti 



7 



7 
'48m 



128m 



(n ■ I/) [i/ ■ (n X k)] (i/ ■ k) QTiT^, 



128m 



^[(i/xn)-p](k-n)f7r3^r,, 

^[(,.xn).k](k.n)nr3^r,. 

^[(:.xn).(i/xa)](k.n)r3V.. 



3fifc2 

7 



m^^^^'y^^ 



Qlkg.^lT 



,2^4 



-[(k.a)-(k.i.)(a.i.)]n^T3* ^ 

"^ ;,n2n2 

(l/ X k) • pn^Tg" 



--kO^nlRTi 



16m 16 8 



'32m^ " 32m^ "" 16m 



7 4 -g^feg^rsr^^ T^%^r3^r.. 

— -k-aTin, —^ l^^ 

-|- (n ■ I/) (k ■ n) (i/ ■ p) r|r,, A^2fc„^r|r,, 

Ibm, Id 



|^(n..)(k.n)(..k)T|T. £.^T|T. 

^ (n ■ z.) (k ■ n) {u ■ a) T*T.. ^^Ikg.TiT.^ 

- -^ (n ■ I/) [i/ ■ (n X k)] (i/ ■ p) OTa'T,, 
b4m 



^ [(i/ X n) ■ (i/ X p)] (k ■ n) T|r.. ±j.2^^,^y3y^^ 

^[(z.xn).(:.xk)](k.n)T|T.. 

-^ [(z. X n) ■ (z/ X p)] [1/ ■ (n X k)] nTfT,, 
'^^ [(i/ X n) ■ (z/ X k)] [i/ ■ (n X k)] nra'T;,^ 



3;ifc2 

32m 


-l- 


Tin^ 






























kiy 

128 


■u 


.Tin, 


128'''^'' 


n' 


RTin, 
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For other values of pulses' areas, using Eqs. H48() one arrives at the expression for the excitation probability 
consisting of 10 terms. Among them five terms associated with transferring coherences between different pairs of 
pulses correspond to Ramsey fringes |22| and play no roll on Doppler broadened transitions |23| ■ Another four terms 
are originated from the population transfer between fields '24|; they comprise the background of the signal. Only one 
term is responsible for the atom interference. Piecing together all these terms one arrives at the following expression 
for the excitation probability: 



1 



2 ^2 

1 — COS di cos 6*2 COS 6*3 — sin 6*3 sin — si 



O2 

— sin6li COS (0„ -(j)ji-(p2-(j)i + 2(j)2) 



where 9j is an area of the pulse j. 



VII. FOUR-PULSE ATOMIC SENSORS 



(57) 



Consider now the interaction of atoms with the sequence of four pulses 

(7r/2, Ti = 0, ki) - (tt, T2, k2) - [tt, T^, k3) - (7r/2, T^, \^i) . 
Using Eqs. (|21|) . H49|) . and H5U|) one finds the upper state distribution 

p,, (r, p, T4 + r) = i {/ [R (R (R (r, p, T3 - r4) , P (r, p, T3 - T^) - ?ik3, T^ - n) , 

P (R (r, p, T3 - T4) , P (r, p, T3 - T4) - hk3, T2 - T3) + nki, -T2) , 
P (R (R (r, p, T3 - 7^4) , P (r, p, T3 - T^) ~ ^kg, T^ - n) , 
P (R (r, p, T3 - T4) , P (r, p, T3 - Ti) - Tikg, T2 - T3) + ^ki, -T2) - ?iki] 
+ / [R (R (R (r , p - ?ik4 , T3 - T4) , P (r , p - ?ik4 , T3 - r4) + ?ik3 , Ts - T3 ) , 
P (R (r, p - ?ik4, n - T4) , P (r, p - ;ik4, T3 - Ti) + hk3,T2 - T3) - ?ik2, -T2) , 
P (R (R (r, p - ?ik4, n - T4) , P (r, p - ?ik4, T3 - T4) + Tikg, Ta - T3) , 
P (R (r, p - ?ik4, n -Ti),P (r, p - ?ik4, T3 - T4) + Tikg, T2 - T3) - ?ik2, -T2)]} 

+ I cos [k4 • r - 2k3 • R (r, p, T3 - T4) + 2k2 • R (r, p, T2 - T4) - ki • R (r, p, -T^) 



-di2 {Ti - 2T3 + 2T2) - ((/.4 - 2(/)3 + 2(/.2 - (/-i)] / 



R(r,p,-T4),P(r,p,-r4)-^ 



(58) 



(59) 



Performing integration over the phase space (using arguments of the / function as integration variables in each term) 
one arrives at the following excitation probability: 



^Wbg+w drdp cos (</>„ -(I>r-(I>4 + ^cj)^ - 2(/«2 + (t>{) f (r, p) , 



Wbg=W= -, 



U--R(r,p + ^,T4] ^2k3-R(r,p+^,Ts) + 2k2 • R ("r, p + ^, Ts 

,5i2 (T4 - 2r3 + 2T2) , 



(60a) 

(60b) 

ki • r, (60c) 

(60d) 



where Whg and w are the background and the amplitude of the interference term correspondingly. Interferometer 
phase H60cll coincides with that given in Eq. |(5J if all pulses have the same effective wave vectors, k^ = k. 
For arbitrary areas of pulses one gets 34 terms, among them are the following: 

• Eight terms correspond to the transfer of level populations between fields, they comprise background of the 
excitation probability given by 



Wbg = - (1 — cos 6*1 cos 62 cos 6*3 cos 64) 



(61) 



• 17 Ramsey terms, which are washed out after averaging over momenta. 
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• Six terms correspond to the coherence transfer between two adjacent time intervals. They have structure similar 
to the case of the excitation by three pulses and are also washed out at least because the ratio of time intervals 
does not satisfy echo conditions. 

• Two terms w± associated with the stimulated echo. They are given by 

w± = — .25sin6'4sin6'3 sin6'2 sin6'i / drdp 

r Iff hki \ f hki \ hk2 

xcos|k4-R RfRfr,p+^,r2J ,P(r,p+^,T2J ±^,T3-T2 
P(R(r,p+^,r2],Pfr,p+^,T2J±^,T3-r2]T^,T4-r3 



-k, R 



R[r,p + ^,T2J,Pir,p+^,T2J+^,n-T2 



5i2 {Ta - T3 



hki 



k2 -R ( r,p+ —^■,T2 ] -ki -v -512T2 - v^j t" v^i 



/(r,p); 



(62) 



A term associated with the double- loop interferometer. It coincides with the second term in Eq. (|60a|l . except 
that the amplitude is equal to 



1 2 L'3 2 ^2 

w — — sin 6i sin — - sin — - sin 9i . 

2 2 2 



(63) 



A. Atomic gyroscope 



This is evident from Eq. (|60c() that, for an arbitrary rotation and acceleration, one can eliminate phase shifts of 
the order T" only if wave vectors and ratios of time intervals (|3J) satisfy the equality 



k4 - 2k3t?j' + 2k2t^ - ki^l = 0. 



(64) 



We always need to eliminate zero-order terms [phase matching condition Hla|) ] and first-order terms (to get interfer- 
ence). For the gyroscope one needs also to eliminate the second-order term. From the system of Eqs. H64|) with n = 0, 
1, and 2 one concludes that all effective wave vectors have to be collinear. 



kj — Sjk, 



(65) 



where k is comprised of counterpropagating traveling waves, Sj = sina^, and aj is one-half of the angle between 
optical fields wave vectors qij and q2j presented in the pulse j. Resolving the set of equations, for example, with 
respect to variables si, S2, and S3, one finds that 



Sl = S4 



S2 



S3 



(1 


~ t2) (1 - 


' h) 


Si 


t2h 
1- h 




2 

Si 


t2{h- 

1- t2 


hV 



2 t3 ( ^3 - t2) 



(66a) 
(66b) 
(66c) 



One sees that for given time delays between pulses one can always choose proper values of effective wave vectors to 
construct a double-loop interferometer. Moreover, wave vectors can be scaled by a given factor. Since the interferom- 
eter phase ()60c(l is linear with respect to wave vectors and, therefore, proportional to this scale factor, one can choose 
it from conditions 



s, < 1, 



(67) 



I.e., 



Si = mm 



i2 t3 



2t2{t3- is) 2 t3 ( t3 - t2) 



(1- t2)(l- h)-^ 



h 



t2 



(68) 
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This expression guaranteed that all wave vectors satisfy the condition H67() and at least one of them yields the maximum 
value k (or Sj = 1). 

To get phase shifts of the atomic gyroscope, one can simply notice from comparison of Eqs. I)54b|l and H60c() that 
it is sufficient to replace T^ in column 1 of Table Q] by 



^3 ~^ ^A 



Si - 2s3^g + 252^2 

1 -21-" ■ 



(69) 



Since all T terms are eliminated, for the launch momentum 

p <C rna/n 

the main phase shift is presented in the fourth row of Table Q] Saving only this term after the replacement H69|) one 
finds 



V 



f2 • (k X a) Ti^P, 
'-±{1- i2)(l- h). 



(70a) 
(70b) 



Phase dependence (I70|l on the relative position of the second and the third pulses (on the parameters ^2 and t^) for 
the optimal value (|68|l of the parameter S4 is plotted in Fig. [3 



[kxal-OT/ 




(T3-T2)/(T4-T2) 



FIG. 2: Dependence of the gyroscope phase on the position of the second and third pulse. 



One can see that the phase (|70f) has only one maximum and that this maximum occurs for the double-loop inter- 
ferometer shown in Fig. ^a), when 



1 and t2 = 1/4, ig = 3/4, 



(71) 
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and (/)^ is given by Eq. 10. 

Nevertheless, one has to move out of the point H71|l at least slightly to wash out contributions from the stimulated 
echo Hfi2|l . Let us consider this case. If 

s, = l- a, (72) 

and aj <C 1, then from Eqs. (|66|l one finds with accuracy up to second order terms in CTj, that the second and the 
third pulses have to be located at 

1/^9 13 3 1 9 2 2 3 2\ ,r,n ■. 

i2 ~ I I 1 + gO-i - era - -(74 + |cricr2 - Yg<^i'^4 + |0-2cr4 + ^cTi - (^2 ~ ^'^tj > (73a) 

1 /o 3 5 1 1 5 15 2 2 13 2\ 

ts ~ ^ ( 3 - -<Ti + (72 - gcr4 + -(72(74 " Jg'^^'^'^ ~ 1^^'^'^ ^ 32^^^ + '^2 " ^(^4] ■ (73b) 

We found these values from Eqs. (|66a|l and (j66b|l . Substituting them into Eq. ()66cf) one gets the phase matching 
condition, in the first order over aj , 

(74 - 2(73 + 2(72 - (71 = 0, (74) 

and the equality in the second order. We calculate, with the same accuracy, the main term (|70ll in the gyroscope 
phase, 

11 1 45 2 1 2 13 2 3 3 1 

W ~ (74 ay a^ (7o CT^ H (7lCT2 (71(74 ^ (79(74, 75 

^ 16 48 24 ^ 1024 ^ 16 ^ 1024 ^ 32 ^ ^ 512 ^ ^ 32 ^ ^' ^ ' 

while all terms in the zero-order over aj, calculated using Tableland the rule (|69|l . are presented in Table UTI In the 
last two columns of this table and of Table HTll below we calculated phase for Cs interferometer (m = 2.21 x 10~^^ g, 
k = 1.47 X 10^ cm-i) for g^ = -980 cm/s^, T4 = 0.8 s, R= 6.72 x 10^ cm, latitude 41°. 

B. Atomic accelerometer 

One can use the four-pulse interferometer to eliminate third-order contributions to the interferometer phase in 
order to increase the accuracy of acceleration measurements. In this case wave vectors and relative positions of the 
second and the third pulses have to satisfy Eqs. (|64() for n — 0, 1, and 3. Resolving these equations with respect 
to wave vectors of the first, second, and third pulses and requiring that at least one of the pulse is comprised of 
counterpropagating waves, one finds 

(1-^3) (1-^2) (1+^2 +^3) ,„. . 

S2 = Trm2 — ^' (76b) 



*3 = -TTTTa HT' (76c) 



S4 = mm 



2 ^3 (^3 — ^2) 

^2^3(^2 + ^3) 2t2{tl-tl) 2h{e^--tl) 



{i-h){i-t2){i + t2 + hy i~tl ' i-tl 



(76d) 



To get the main contribution to the interferometer phase one should perform the replacement (|69|l for n = 2 in the 
very first cell of Table P It yields the following expression: 

2{t2+t3) ^ ^ ' 

This dependence is shown in Fig. 

One sees that only one maximum of the accelerometer 's phase corresponds to values of ^2 and ^3 given by Eq. I^. 
Substituting these values into Eqs. (|76l) one verifies that the maximum of accelerometer phase occurs when all pulses 
are comprised of counterpropagating fields, Sj = 1. 

From Eq. @ one can see that T4 — T3 ^ T2 and therefore the contribution from stimulated echo processes given 
by Eq. (|62|l is washed out after averaging over momenta. 

Different terms in the accelerometer phase can be obtained from Table ^ using the replacement (|69|) . They are 
presented in Table ITTll 



16 



TABLE 11; Atomic gyroscope phase shifts; columnl, the general case; column 2, for the wave vector k and launch momentum 
p in the west-east and vertical directions; columns 3 and 4 represent the absolute and the relative phase shifts when p = for 
Cs gyroscope 
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VIII. DISCUSSION 



In this paper we showed that, in the presence of inertial forces and homogeneous acceleration gradient, the density 
matrix in the Wigner representation in the free space still obeys classical Liouville equation for the distribution 
function. It allows us to get interferometers' phases without calculating path integrals. 

For all interferometers under consideration we calculate phase shifts up to terms proportional to T^ for an arbitrary 
constant in time acceleration and rotation. 

Using double-loop interferometers 7r/2 — tt— tt — 7r/2 one can eliminate the phase dependence on the acceleration only 
and get an interferometer with the main phase term proportional to the cross product of acceleration and rotation 
frequency. One can use this interferometer as a gyroscope to measure rotation frequency components perpendicular 
to the acceleration. Even for the rotation frequency as small as the Earth rotation rate, the phase of this gyroscope 
can achieve (f>m ~ 1^)0. Eliminating the acceleration term, which is typically five orders of magnitude larger, could 
be incomplete if pulse areas arc not precisely equal to tt/2 or tt, because in this case stimulated echo processes, 
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0.06 



(k-a)T4^ 




(T3-T,)/(T4-T2) 



FIG. 3: The same as in Fig. |5| but for the phase of accelerometer. 



which are extremely sensitive to the acceleration, contribute to the probability of the particles' excitation. To exclude 
these processes we propose to drive Raman transitions between hyperfine atomic sublevels by a pulse of travelling 
noncounterpropagating optical waves. For this goal one can use even slightly noncounterpropagating fields having 
effective wave vectors given by Eqs. ()f)5|l and H72|l . where the small parameter aj determines the deviation of the pulse 
j wave vector from the maximum value k. If wave vectors satisfy the phase matching condition Hla|) for double- loop 
interferometer but violate the phase-matching condition ljlb|l for stimulated echo, i.e., if 



^dl ^ (7i — 2(73 + 2(72 - CTl = 0, 
(^st = 0-4 - (T3 - 0-2 + O"! 7^ 0, 



(78a) 
(78b) 



integrands in Eq. H62I) for probabilities of excitation associated with stimulated echo processes become rapidly os- 
cillating in space function, which have the period of the order of A/ust, where A = 2tt /k is the effective wavelength. 
When 



A 



^cloud 



< (Tst < 1, 



(79) 



where Scioud is the initial atomic cloud size, in the case when the initial atom density is smooth in space, probabilities 
(|62|l become exponentially small with respect to the parameter UstS cloud /^ ^ 1. In addition to arguments related to 
the difference in phase matching conditions, one can expect further decrease of probabilities (|62|l due to the fact that 
pulses are not properly positioned in time [1 — ^3 7^ t2, where tj are given by Eqs. H73() ] to produce the stimulated 
echo. We obtained the required pulses' positions and the main term in the gyroscope phase for arbitrary ctj, see Eqs. 
(|73|l . (|7()a|l . and (|75|l . But it is sufficient for only two effective wave vectors to have ctj 7^ 0. If, for example. 



C3 = (74 = and CT2 



0-1 



^0, 



(80) 
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TABLE III: The same as in TableQ] but for the four-pulse interferometer shown in Fig. H^c). The last two columns are for the 
phase shifts' numerical and relative values of the Cs interferometer for p = 0. 
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8-3^5 2 „ ^4^ 
^-^ \{v X n) . (i. X a)] (k . n) TtT^^ "— = 

128 ^ 

^ 

then a St — — , second and third pulses have to be apphed at moments 

^4 /^, , 5 , 13 2 



and gyroscope phase is given by 












T.^^\l+ 3^1 + ^.f ) , (81a) 

Ta^rJl + l.^ + l.?) (81b) 



A "side-effect" of using Raman pulses with nonequal effective wave vectors is that in this case one gets Ramsey 
fringes, i.e., excitation probability oscillating dependences on the Raman detuning 5i2 having period of the order 
of the inverse time delay between pulses T~^ [see Ramsey phases 0^ (I54c|) and (|60d|l for single- and double-loop 
interferometers] . Evidently, this new type of Ramsey fringes arises because even for Raman pulses with copropagating 
effective wave vectors, in the condition of the Doppler phase cancellation, time separations between pulses Tj enter 
with weight factors kj , while these factors are absent in expressions for the Ramsey phase. In contrast to the previously 
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proposed ^3| and observed ^ Raman-Ramsey fringes on Doppler broadened two-quantum transitions based on the 
use of Raman pulses having counterpropagating effective wave vectors, Ramsey fringes found here do not undergo 
recoil splitting, i.e., they are centered at the two-quantum line center. 

Double-loop interferometer shown in Fig. ^c) allows one to get the phase where all cubic in T terms are eliminated 
and can serve as an accelerometer. We expect that owing to this cancellation, the accuracy of the acceleration 
measurement can be increased. 

Thus in this paper we propose the following scenario for the atomic gravity and rotation sensing. First one measures 
acceleration a using three mutually perpendicular double-loop accelerometers. Afterwards, one uses the spatial-domain 
atomic gyroscope Q , orients it in the plane perpendicular to a and measures the component of the rotation frequency 
ri along a. Then one uses two mutually perpendicular time-domain double-loop gyroscopes to get components of J7 
perpendicular to a. 

Further obvious step here could be the consideration of multiple- loop interferometers. In general for an interfer- 
ometer consisting of (, unidirectional pulses (first and last are 7r/2 pulses, others are tt pulses), having the same wave 
vector, one has £ — 2 parameters tj = Tj/Ti (j = 2, . . . , ^ — 1), which can be chosen as a root of the system of ^ — 2 
equations 

1 - 2t;i^ + 2t;^2 + • • • + (-1)^ 2^2^ = (83) 

to eliminate phase terms evolving as T"^, . . . ,Tg '^^ . In the array u = (mi, . . . , M£_-2) one of parameters Uj has to be 
equal to 1, to eliminate the Doppler phase, while the choice of other parameters depends on phenomena one would 
like to sense using the interferometer. Remaining noneliminated phase terms evolving as T^ can be obtained from 
Table m by replacement of 

n '-^ ^^^^ ^ ^ ' t/. (84) 

The case of four-pulse interferometers is considered above in detail. Using five-pulse interferometers one can eliminate 
T^ dependencies. For gyroscope u = (1, 2, 4) , we found the numeric solution of system (|83|l . 

i2 = 0.16615, ^3 = 0.54057, U = 0.87442. (85) 

This gyroscope phase is given by 

0„ = -1.0143 X 10^2 j^ • (k X a) r| + O (T|) . (86) 

For the Cs gyroscope with the same parameters as in Table HTl 0^ ~ —41.2. 
For the five-pulse accelerometer, u = (1, 3, 4) , one finds 

t2 = 0.20594, tg = 0.594 37, U = 0.88843, (87) 

and the phase is given by 

0„ = 2.1557 X lO^^k • aTl -I- O {tI) . (88) 



For the Cs accelerometer with the same parameters as in Table UTTl 0„ « —2.00 x 10^. 

Another example of five-pulse interferometer is a "figure 8 1/2" interferometer proposed in [lj| to eliminate T^ and 
Tg terms, such as the dominate contribution to the phase is caused by the infiuence of space-time curvature. 

Consider also six-pulse interferometer. One can use this pulse sequence to eliminate all phase terms arising as a 
result of the atom motion to the accuracy Tg choosing u = (1, 2, 3, 4) . The numerical solution of the system H83I) is 
given by 

i2 = 9.5492 X 10-^ H = 0.34549, U = 0.65451, tg = 0.90451. (89) 

Of course, one cannot use this interferometer to sense atomic motion. But one can use it for the precise measurement 
of the other spectroscopic data. The accuracy of this measurement improves when the time of evolution Tg increases. 
Earth gravity and rotation superimpose limitation on this time. This time can be increased in the microgravity envi- 
ronment. We actually propose to use the six-pulse interferometer instead of transferring experiments into microgravity 
environment. 

Possible examples here are measurement of levels polarizability 26] or observation of the Aharonov-Bohm effect 
|27j . We also believe that one can insert additional tt pulses in the scheme of the recoil frequency measurement 6] to 
make this measurement independent on the gravity, gravity gradient and the Earth rotation. 
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(a) 



T2 '' T3 T4 Ts 

7i/2 n n n nl2 



(b) 



T2 T3 T4 T5 

7l/2 71 71 71 7l/2 



(C) 




7l/2 71 71 71 71 7l/2 



FIG. 4: Recoil diagrams for the five-pulse gyroscope (a), accelerometer (b) and six-pulse interferometer insensitive to both 
rotation and acceleration. 

Recoil diagrams for five- and six-pulse interferometers are shown in Fig. 01 Regarding the six-pulse interferometer, 
one sees that it has to be protected from two kinds of the stimulated echo processes produced by the sequences of 
pulses {1, 2, 5, 6} and {1, 3, 4, 6} . A similar note can be made regarding the "figure 8 1/2" interferometer 14], where 
both stimulated echo produced by pulses {1, 2, 4, 5} and one-loop interferometer produced by pulses {1, 3, 5} have to 
be excluded. An effective technique of eliminating all these unwanted signals is using Raman pulses with noncqual 
effective wave vectors, as we demonstrated and analyzed above for the double-loop atomic gyroscope. 
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